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Abstract 

Between 2008/3/24 and 2008/4/2, the three active regions NOAA ARs 10987, 
10988 and 10989 were observed daily by the Synoptic Optical Long-term In- 
vestigations of the Sun (SOLIS) Vector Spectro-Magnetograph (VSM) while 
they traversed the solar disk. We use these measurements and the nonlinear 
force-free magnetic field code XTRAPOL to reconstruct the coronal magnetic 
field for each active region and compare model field lines with images from the 
Solar Terrestrial RElations Observatory (STEREO) and Hinode X-ray Telescope 
(XRT) telescopes. Synoptic maps made from continuous, round-the-clock Global 
Oscillations Network Group (GONG) magnetograms provide information on the 
global photospheric field and potential-field source-surface models based on these 
maps describe the global coronal field during the Whole Heliospheric Interval 
(WHI) and its neighboring rotations. Features of the modeled global field, such as 
the coronal holes and streamer belt locations, are discussed in comparison with 
extreme ultra-violet and coronagraph observations from STEREO. The global 
field is found to be far from a minimum, dipolar state. From the nonlinear models 
we compute physical quantities for the active regions such as the photospheric 
magnetic and electric current fiuxes, the free magnetic energy and the relative he- 
licity for each region each day where observations permit. The inter connectivity 
of the three regions is addressed in the context of the potential-field source- 
surface model. Using local and global quantities derived from the models, we 
briefiy discuss the different observed activity levels of the regions. 



1. Introduction 

A large part of the low solar corona is dominated by the magnetic field. This field 
is created inside the Sun and then emerges into the atmosphere. The magnetic 
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field plays a defining role in most coronal phenomena, quasi-static structures 
from cool, dense prominences to hot loops, and dynamical events such as flares 
and coronal mass ejections (CMEs). The most spectacular events with the great- 
est impact on us originate from the lower corona. If we are to understand the 
various phenomena that occurred in the corona during WHI it is therefore crucial 
that we determine the 3D coronal magnetic fleld as well as possible. Unfortu- 
nately, the coronal fleld is notoriously difficult to measure ( [Tomczyk et a/., 2007 



|Liu and Lin, 2008] ), signiflcantly more difficult than in the much denser pho- 
tosphere. It has therefore become common practice to reconstruct the coro- 
nal fleld from boundary conditions derived from photospheric fleld measure- 
ments ( |Amari et a/., 1997| [Schrijver et a/., 2006[ [Wiegelmann, 2008] ). Our best 



estimate of coronal magnetic structure comes from such modeling efforts. 

On the other hand, it is clear that there remain significant problems with the 
accuracy of the models. The governing magnetohydrodynamics (MHD) equa- 
tions are complex and even their equilibrium solutions are difficult to constrain. 
Because the magnetic field so dominates the corona, most coronal fields can 
reasonably be assumed to be approximately force-free. Models of the global 
lower coronal field usually do not include electric currents in the lower coronal 
domain either. Such models, referred to as potential-field source-surface (PFSS) 
models, incorporate on their upper boundary a current "source surface" where 
the magnetic field is forced to be radial, mimicking the effect of the solar wind 
dAltschuler and Newkirk, 1969||Schatten, Wilcox, and Ness, 1969||Hoeksema, 1984H 



|Wang and Sheeley, 1992D . The global field produced by such models usually 



closely approximates much more expensive MHD models in practice ( [Neugebauer et a/., 1998( 1 



|Riley et al.^ 2006) and they have been successfully applied in studying basic 
properties of the global coronal field such as coronal holes, streamer belt locations 
and dipole axis tilts (Zhao, Hoeksema, and Scherrer, 2005} [Schrijver and Liu, 2008|).| 



Although [Schrijver et al (j2005p found that many quiescent active regions 
resemble potential field models, recent attempts to measure the accuracy of 
coronal magnetic field models extrapolated from photospheric field measure- 
ments have produced discouraging conclusions. Since the launch of NASA's twin 
Solar Terrestrial Relations Observatory (STEREO) satellites it has been possible 
to estimate from simultaneous pairs of extreme ultra-violet (EUV) images from 
distinct vantage points the path through 3D space of any loop structure that can 
be correctly identified in both images. fSandman et a/.l(|2009p used this technique 
to test the accuracy of potential-field models of three active regions from line- 
of-sight MDI field measurements. They measured the accuracy by calculating 
the misalignment angle between each observed loop and the model at different 
positions along the loop and noting the median value. For the three regions 
studied they found averages and standard deviations 25° ± 8°, 19° ± 8° and 
36° ± 13°. We remark that these observational tests do not involve direct com- 
parison between observed and modeled coronal fields. It is not known how closely 
the loop structures observed in EUV adhere to magnetic field trajectories. 

As for NLFFF modeling, the picture is still more complicated. There are 
several very different solution methods available and, unlike the potential-field 
problem, the nonlinear force-free problem cannot be immediately cast in a clear, 
mathematically well-posed way. The force-free equations are nonlinear and of 
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mixed mathematical (elliptic/hyperbolic) type. Only two pieces of information, 
e.g., the normal components of the field B and the electric current J, can be 
imposed on the boundary to give a mathematically well-posed boundary value 
problem, whereas observations of all three magnetic field vector components are 
available. Different solution methods apply this information in different ways, 
often producing confiicting results. To add to the complication, Grad- Rubin 
methods ( |Grad and Rubin, 1958D such as the one used here can incorporate 
data from either the positive or the negative magnetic polarity but not both in 
a single model. Models from the two polarities can be very dissimilar, and often 
are in practice. On the other hand, sometimes it is clear which polarity one must 
apply as boundary data ( |Canou and Amari, 2Q1QD . New approaches to take into 
account both polarities have been recently proposed by [Wheatland and Regnier 
|2Q Q9^ and [Amari and Aly| ([20101). 

dSchrijver et a/., 2008D and ( [DeRosa et a/., 2009D compared several NLFFF 
codes using common boundary data sets comprising Hinode vector magnetograms| 
embedded within line-of-sight MDI magnetograms. They found major disparities 
between the various models, in their magnetic field and electric current struc- 
tures. dSchrijver et a/., 2008 ) find that the model fields differed significantly in 
geometry, energy content, and force- freeness. ( DeRosa et a/., 2009[ ) concluded 
that while models can coincide well with observed features it remains difficult 
to determine whether a significant proportion of a modeled field accurately re- 
produces the coronal field over the entire coronal volume above an active region. 
dCanou and Amari, 2010D have re-analyzed the data set studied by DeRosa et 
al. (2009), modeling the field of AR 10953 exploiting the fact that their Grad- 
Rubin method (also used here) takes its information on the photospheric electric 
current (via the force-free function a) from only one of the magnetic polari- 
ties. They chose to integrate from the polarity most fully observed by Hinode, 
including a intense sunspot electric current fiux, thereby avoiding the lack of 
information on the electric current in the MDI portion of the boundary data set. 
Finding evidence of a bald patch in the vector magnetogram they successfully 
modeled a magnetic fiux rope that corresponded well with a Hinode XRT sigmoid 
and a filament seen in SMART Ha data. 

A further complication is that photospheric vector field measurements of 
a given solar phenomenon show significant disagreement from instrument to 
instrument and even between different inversions of a single Stokes data sel0. 
In view of this, it may be optimistic to expect a coronal model extrapolated 
from photospheric data to agree with coronal observations. Wiegelmann et aL\ 
d2010p show encouraging evidence that nonlinear force-free field models based 
on good photospheric vector field measurements can give reasonable estimates 
of the structure of the coronal field. 

Nonlinear force-free models give estimates of the free magnetic energy and 
magnetic hehcity of active region structures. Such modeling techniques can 
capture twisted coronal magnetic structures, providing evidence that twisted 
fiux ropes exist in equilibrium in the solar corona (e.g.. ICanou "eTaT] ([20091) and 
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sometimes detailed knowledge of the different twisted structures in an active 
region can clarify discrepancies between the handedness of twisted flux tubes 
ejected from the Sun (magnetic clouds) and the handedness of twisted flux low 
in the corona dRegnier and Amari, 2QQ4t . Usually electric currents flow in more 
than one direction in an active region. Nonlinear force-free fields can capture 
this current structure and allow a more accurate and detailed estimate of active 
region helicities. 

In this paper we combine the two different types of coronal model discussed 
above, PFSS and NLFFF, to describe the global coronal field structure and the 
active region field properties during the WHI. The PFSS models will be based 
on Global Oscillations Network Group (GONG) synoptic magnetograms and the 
NLFFF models extrapolated from Synoptic Optical Long-term Investigations of 
the Sun (SOLIS) vector magnetograms. From the PFSS models for the WHI 
and its neighboring rotations we can study the coronal hole distributions and 
streamer belt locations, comparing them to synoptic EUV and coronagraph 
data from NASA's STEREO spacecraft. Meanwhile NLFFF models of the active 
regions can provide information on the magnetic field including its basic struc- 
ture, free magnetic energy and relative magnetic helicity. We will compare the 
modeled field structures to Hinode/XRT and STEREO/SECCHI/EUVI images 
and interpret the PFSS and NLFFF models together to characterize the active 
regions. 

The paper is organized as follows. We will describe the modeling techniques 
in Section [2j Then, we will present PFSS models for Carrington rotations 2067-9 
and associated observations of the global coronal structure in Section [3l before 
presenting NLFFF models and observations of the active region structure in 
Section |4] including a brief study of the interconnectedness of the regions and 
their temporal evolution. We will conclude in Section [5l 

2. The PFSS and NLFFF models 

In the low corona, typical Alfven speeds are so fast (1000 km/s) compared to 
the observed evolution time of structures that any persistent structure is static 
to a good approximation. The magnetic field dominates the dynamics such that 
the Lorentz force J x B, unopposed by significant plasma forces, is negligible 
such that almost all Maxwell stresses are contained within the field and the field 
and electric current are approximately parallel, 

J = a(x)B. (1) 

Because the quantity ol is inversely proportional to the length scale, the current- 
free approximation a ^ is useful for studying large-scale structure. 

Under the simplifying assumption that there are no significant electric cur- 
rents in the corona, the coronal magnetic field may be represented by B = 
— V?/^(r, ^, (\)) in spherical coordinates where ?/^(r, ^, (/)) is a spherical harmonic se- 
ries in d and (\). The spherical harmonic coefficients are chosen such that the radial 
field component at the lower boundary of the model matches the photospheric 
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magnetogram while the r-dependence of ijj guarantees that the field is radial at a 
"source-surface", the outer surface of the model ( [Altschuler and Newkirk, 1969 ; 
[Hoeksema, 198 4^). While line-of-sight boundary data can also be used, radial 
boundary conditions derived from line-of-sight data generally give better models 
because the photospheric field is approximately radial and because polar fields 
are poorly represented by line-of-sight measurements ( |Wang and Sheeley, 1992[ 
|Petrie and Patrikeeva, 2QQ9[ ). For a given magnetogram this method in theory 
provides a unique solution that is guaranteed to exist. In practice the solution 
depends not only on the magnetogram, but also on the location and shape of 
the source surface and the number of spherical harmonics used. For the PFSS 
models presented in this paper we use 36 spherical harmonics and we apply a 
spherical source surface of radius 2.5 solar radii. The GONG synoptic maps used 
have 360 x 180 pixels in longitude-sine (latitude) coordinates. 
A force-free field obeys the equations 



V X B = aB, (2) 
V-B = 0. (3) 



The function a(x) must satisfy 



B . Va = 0. (4) 

Equations (j2l H]) have mixed elliptic- hyperbolic structure, consisting of an 
elliptic problem for B with fixed a and a hyperbolic problem for a with fixed 
B. The nonlinear force- free field calculation is of Grad- Rubin type, i.e. B and 
a are calculated iteratively with the elliptic and hyperbolic parts of the prob- 
lem solved successively at each step. To solve the elliptic problem one needs 
to impose the vertical field component on the boundary, whereas solving 
the hyperbolic problem requires information about a on either the part of the 
boundary where 5^ > or the part where 5^ < 0. Although all three com- 
ponents of the photospheric field are measured, only two pieces of information, 
e.g., the vertical components of the magnetic field and the electric current, may 
be imposed on the boundary to obtain a mathematically well-posed boundary- 
value problem for the mixed elliptic- hyperbolic NLFFF equations. For this type 
of well-posed formulation, rigorous existence and partial uniqueness theorems 
have been proved ( [Bineau, 1972[ |Boulmezaoud and Amari, 2000D . In addition, 
the finite- difference XTRAPOL code used here works with a vector potential A 
associated with the magnetic field, such that B = V x A, and with a staggered 
mesh such that the approximate solution is a member of the kernel of the V- 
so that the model field is constrained to be divergence-free to machine precision 
( [Amari, Boulmezaoud, and Mikic, 1999[ ). More recently, the code was extended 
to accommodate solutions with non-zero B and a on all six sides of the compu- 
tational box ( [Amari, Boulmezaoud, and Aly, 20061 ). ^^i^ paper we exploit all 
of these properties of the code. All of the NLFF calculations in this paper were 
performed on 120 x 120 x 80 grid points using a non-uniform grid where grid 
points are accumulated in region of strong current and strong magnetic field. 
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In comparison a typical resolution of the SOLIS AR vector magnetograms is 
300 X 300 grid points. For each magnetogram we calculate a model from each 
polarity. If there is significant disagreement between them there is usually a sim- 
ple explanation, as we discuss in later sections. The NLFF calculation requires 
a potential field solution. This is calculated independently of the global PFSS 
model discussed above, in the same Cartesian domain as the NLFFF model, 
using the SOLIS measurements for the vertical field component as boundary 
conditions. The instrumental noise level of the VSM is a few G for the longi- 
tudinal field and about 40 G in the transverse direction (Jack harvey, private 
communication). The results of the Stokes inversion calculation are not reliable 
for fields weaker than about 50 G. In practice we include only fields stronger than 
about 50 G in our boundary data. Because of the enormous differences in spatial 
resolution between the GONG synoptic map and the SOLIS magnetograms for 
the active regions, the question of consistency between these data sets and the 
associated coronal field models is not a trivial issue. We will see as the paper 
progresses that, while the gross features of the active regions appear consistently 
in both data sets, the distinctive characteristics of the three regions are much 
more obvious in the vector data. 

3. The global structure 

Figure [1] shows the GONG synoptic maps for Carrington rotations 2067-9. These 
show that rotation 2068 had three well-defined active regions, NOAA Active 
Regions 10987, 10988 and 10989 at approximately equal latitude, about —10°, 
and separated by about 30° of longitude from each other; 10989 at about 205°, 
10988 at 235° and 10987 at 265°. Ah are consistent with Hale's polarity law in 
having leading negative polarity but their tilt angles appear to differ. The maps 
also suggest that AR10989 was shorter-lived than the other two regions. It is 
dispersed in the 2069 map while another less intense region (not numbered by 
NOAA) appears at the same Carrington longitude but North of the equator and 
with opposite leading polarity. Meanwhile regions 10987 and 10988 reappear in 
rotation 2069, both showing evidence of wear and tear from turbulent diffusion. 
The map of rotation 2067, on the other hand, shows little coherent structure but 
only diffuse field around longitude 235°, latitude —10° where 10988 will appear. 

Figure [2] show features of the PFSS model for rotations 2067-9 based on 
the synoptic map of Figure [H These plots show the main neutral line which is 
formed by the apexes of the set of tallest closed field lines projected radially 
downwards onto the photosphere. The coronal streamer belt is expected to be 
located immediately above this neutral line, separating open field regions of 
opposite polarity. The thin black lines in the top left plot denote locations of 
pseudo-streamers, streamer structures sometimes found above multipolar loop 
systems, separating regions of open field of the same polarity. The principal 
morphological difference between the two types of streamers in coronagraph 
images is that the helmet streamer cusps are located a solar radius or two above 
the solar surface whereas only the long stalks of pseudo-streamers are visible 
above such a height. In the PFSS model the closed portion of the streamer 
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Figure 1. GONG synoptic magnetograms of the radial magnetic field component for 
Carrington rotations 2067 (top), 2068 (WHI, middle) and 2069 (bottom). 
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Figure 3. STEREO/SECCHI/EUVI 171 A synoptic maps for CR 2067 (top pictures), CR 
2068 (middle pictures) and CR 2069 (bottom pictures) from the ahead (left pictures) and 
behind (right pictures) spacecraft. These show active regions as bright regions and coronal 
holes as dark regions. 



must extend to the outer boundary because a current sheet cannot appear 
in the domain of a regular potential field. On the other hand, the cusp of 
a pseudo-streamer can appear at any height in the domain. These are found 



by following Wang, Sheeley, and Rich (j2QQ7p . We define a point on the source 
surface to be a pseudo-streamer if at least two of its neighbors have foot-point 
separation> (7r/12)i?5, where Rg is the solar radius, and are of the same polarity. 
Also shown in the model plots are the predicted coronal holes, regions formed 
by the set of foot points of field lines open to the outer boundary of the model. 
The STEREO/EUVI 171 A maps in Figure [3] show the coronal hole locations as 
dark patches while the STEREO/CORl synoptic maps in Figure H] summarize 
the brightness distribution of the East limb over the rotation. The modeled 
coronal hole and streamer locations match the observations reasonably well. 
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Figure 4. STEREO/SECCHI/CORl synoptic maps for CR 2067 (top pictures), CR 2068 
(middle pictures) and CR 2069 (bottom pictures) from the ahead (left pictures) and behind 
(right pictures) spacecraft. These maps show the streamer brightness at 2.6Rs off the east 
limb, where Rs is the solar radius. The equivalent maps for the West limb appear very similar. 



The three ahgned active regions produce a warp in the neutral hue during 
the WHI. They do this by creating an East- West magnetic moment on a spatial 
scale large enough to influence the global field structure. During much of 2008 
there was activity at approximately these Carrington longitudes, corresponding 
to a major filament system at the same longitudes observed throughout 2008 by 
Patrick Mcintosh - see lWebb et a/.l(|2010p for examples. It is this asymmetric ac- 
tivity pattern combined with relatively weak polar fields ( [Petrie and Patrikeeva, 20Q9| )| 
that results in a global coronal field configuration that is far from dipolar. 
Comparison of the spherical harmonic components of the potential field shows 
why the field is so tilted and non-dipolar. The equatorial dipole component is 
approximately half as strong as the polar dipole, producing a combined dipole 
tilted at about 60° to the rotation axis. Even at the source surface the full dipole 
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Figure 5. The PFSS models of Figure[2]extrapolated from the GONG synoptic magnetograms 
of Figure [H for Carrington rotations 2067 (top left), 2068 (WHI, top right) and 2069 (bottom). 
Positive and negative coronal holes are colored red and green. Streamer-belt fields are plotted 
in blue. 



component is only three times stronger than the quadrupole and octupole com- 
ponents of the field. Figure [5] shows the models for rotations 2067-9 in spherical 
coordinates, viewed from longitude 240°. This viewpoint shows the warp and tilt 
of the streamer belt. The strength and alignment of the active regions during the 
WHI produces a slightly larger warp and tilt in the model for 2068 compared to 

2067 and 2069. However, the displacement of the streamer belt from the equator 
is approximately the same for all three rotations according to Figure [H 

Figure [6] displays in spherical coordinates the photospheric surface field for 
rotation 2068 in greyscale with the streamer-belt field in blue and the positive 
and negative open field in red and green. Figure [71 shows the magnetogram for 

2068 in greyscale with features of the coronal field structure over-plotted. The 
neutral line is plotted in blue and the positive and negative coronal holes are 
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Figure 6. "Hairy ball" plot of the PFSS model of Figure [2] extrapolated from the GONG syn- 
optic magnetogram for OR 2068. The magnetogram is plotted in greyscale. The streamer-belt 
fields are plotted in blue and the open positive and negative fields in red and green. 

represented by regions of red and green dots. The red and green lines represent 
the positive and negative field trajectories that connect the base of the corona to 
where the ecliptic plane intersects the outer boundary of the model. These fields 
are therefore the predicted trajectories of solar wind plasma that will eventually 
reach Earth. Unusually for solar minimum, much of this ecliptic fiux originates 
from midlatitude coronal holes. This is another indication of how far from a 
dipolar state the global coronal field was during this rotation. 

4. Structure of the active regions 

The three main active regions of the WHI have been compared and contrasted 
in terms of their CME productivity and fiare activity. The CMEs and other tran- 
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Figure 7. Plot of the positive (red) and negative (green) ecliptic-plane fields of the PFSS 
model for CR 2068. Also plotted are the streamer-belt neutral line in blue, the positive and 
negative coronal holes represented by red and green dots and the GONG magnetogram in 
greyscale. 



sient phenomena during the WHI have been surveyed by Webb, Gibson, and Thompson! 
(f2QTQ|) , ICremades, Mandrini, and Dasso| (|2QTQ|) and |Sterhng| (|2QTQ|) . These sur- 
veys were carried out independently and used different data sources and analysis 
techniques. I Webb, Gibson, and Thompson (j2QlQp combined these surveys into a 
comprehensive catalog and IWebb et a/.l(|2010p extended this catalog to CR 2067 
and CR 2069. About two thirds of the CMEs from CRs 2067-9 occurred dur- 
ing CR 2068 and most of those between March 21 and April 11. According 
to GOES X-ray flux measurements significantly more flaring activity occurred 
during CR 2 068 than during C R 2067 and CR 2069, mainly between March 24 
and April 6 ( |Webb et a/., 20101 ). This flaring activity was low to moderate. Even 
during the relatively active period of CR 2068 there were only three flares with 
peak fluxes above the X-ray C level (10~^ W/m^): an Ml. 7 flare on March 25, 
by far the most energetic flare during the WHI, and two CI. 2 flares on April 
3. Most of the CMEs were associated with the three ARs 10987-9, in particular 
with AR 10989. As Figure 10 in We bb et~aL\ (|2010p shows, significantly more 
CMEs were traced to locations close to (mainly to the South of) AR 10989 than 
to AR 10987 and AR 10988. The two largest CMEs during the WHI occurred 
on March 25 and April 9. The March 25 event, associated with the Ml. 7 flare 
mentioned above, occurred at the East limb and was associated with AR 10989. 
(The April 9 event was associated with a prominence eruption.) AR 10989 was 
the most CME-productive of the three regions, producing an M2 flare on March 
25th and various ejecta, while AR 10987 was less active and AR 10988 least active 
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Figure 8. SOLIS vector magnetogram (top) and corresponding vertical electric current den- 
sity (bottom) of AR 10987 on 2008 March 27 at 1535 UT. In the plot of the vector magnetogram 
red and blue denote positive and negative vertical magnetic field, with the color bar showing 
the values in Gauss, and the black arrows show the transverse field vectors. In the electric 
current plot red and blue denote positive and negative current with the color bar showing the 
values in statampere/cm-^. 
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Figure 9. SOLIS vector magnetogram (top) and corresponding vertical electric current den- 
sity (bottom) of AR 10988 on 2008 March 29 at 1528 UT. In the plot of the vector magnetogram 
red and blue denote positive and negative vertical magnetic field, with the color bar showing 
the values in Gauss, and the black arrows show the transverse field vectors. In the electric 
current plot red and blue denote positive and negative current with the color bar showing the 
values in statampere/cm-^. 
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Figure 10. SOLIS vector magnetogram (top) and corresponding vertical electric current den- 
sity (bottom) of AR 10989 on 2008 March 31 at 1742 UT. In the plot of the vector magnetogram 
red and blue denote positive and negative vertical magnetic field, with the color bar showing 
the values in Gauss, and the black arrows show the transverse field vectors. In the electric 
current plot red and blue denote positive and negative current with the color bar showing the 
values in statampere/cm-^ . 
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Figure 11. Oblique and top projections of selected representative 3D magnetic field lines of 
the NLFFF model for AR 10987 on March 27. On the lower boundary, the vertical magnetic 
fiux distribution is indicated by red (positive field) and blue (negative field) coloring. The 
colors of the field lines are not significant. 
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Figure 12. Oblique and top projections of selected 3D magnetic field lines of the NLFFF 
model for AR 10988 on March 29. On the lower boundary, the vertical magnetic fiux distri- 
bution is indicated by red (positive field) and blue (negative field) coloring. The colors of the 
field lines are not significant. 
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Figure 13. Oblique and top projections of selected representative 3D magnetic field lines of 
the NLFFF model for AR 10989 on March 31. On the lower boundary, the vertical magnetic 
fiux distribution is indicated by red (positive field) and blue (negative field) coloring. The 
colors of the field lines are not significant. 
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Figure 14. Oblique and top projections of selected 3D magnetic field lines for AR 10987 on 
March 27. On the lower boundary, the vertical magnetic fiux distribution is indicated by red 
(positive field) and blue (negative field) coloring. The colors of the field lines are not significant. 
For this plot the field lines represent the fiux through an electric current concentration in the 
positive magnetic fiux concentration. The resulting field lines together resemble a sigmoidal or 
S-shaped structure. 
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of all. IWelsch and McTiernanI (j2QlQp augmented the standard GOES flare list 
with Sam Freeland's "Latest Events" databas^ and the RHESSI flare catalog 
and found that the activity levels among the three active regions were more 
evenly balanced when the "Latest Events" database was used to associate more 
of the GOES flares to source active regions, and when RHESSI events were 
manually associated with the active regions. 

According to several predictors of flare activity based on the photospheric 
magnetic field (total unsigned magnetic flux, magnetic flux near polarity in- 
version lines, amount of cancelled flux, proxy Poynting flux, relative magnetic 
helicity flux), IWelsch and McTiernanI pUTOl) found that AR 10988 should have 
been the most active region of the three and AR 10989 the least. Of the magnetic 
parameters they studied, the characteristic helicity, the dimensionless quantity 
defined as the cumulative relative magnetic helicity flux normalized by the square 
of the mean magnetic flux over the time the active region was observed, was the 
only one for which AR 10989 had a higher value than the other two regions. 
Magnetic helicity is not easily dissipated in the corona ( [Berger, 1984D and is 
thought to accumulate there until bodily removed by CMEs ( Low, 2001 ). 

We are fortunate that the SOLIS VSM was collecting full-disk vector-field 
measurements daily during the disk passage of the three active regions. We 
therefore have a time series of measurements including an observation close to 
disk-center for each region. While the three regions appear quite similar in the 
synoptic magnetogram for rotation 2068 shown in Figure [U closer examination 
reveals major differences between them. The top pictures in Figures ISlfTOl show 
SOLIS vector magnetograms for the three active regions. Each region is repre- 
sented by the magnetogram closest to central meridian. AR 10987 is a bipolar 
region with approximately equally strong concentrations of positive and negative 
flux. AR 10988 has a strong negative leading polarity, the most intense field in 
the triplet of regions, while the trailing positive flux is significantly weaker and 
more diffuse. Finally, AR 10989 is the weakest and most diffuse region. It has no 
flux concentration comparable to those of the other two regions. These patterns 
are also found in the GONG continuum intensity images (not shown). AR 10987 
has a sunspot of each polarity and they are of approximately equal size. AR 10988 
has a large leading-polarity sunspot, the largest among the three regions, but no 
comparable trailing sunspot. AR 10989 has no significant sunspot. The vertical 
unsigned and net magnetic fluxes, Ffot and F^et, associated with the magne- 
tograms in Figures ISHTOl are recorded in Table [H ARs 10987 and 10988 have 
substantial flux imbalances, over 30% and just under 20%. These imbalances are 
of opposite signs, hinting at some magnetic connection between the two regions. 
Meanwhile, AR 10989 is more nearly fiux-balanced with an imbalance of just 
under 10%. 

The vertical electric current maps in the bottom pictures of Figures ISlfTOl also 
show current concentrations of approximately equal size in AR 10987, a strong 
leading current concentration and weak, diffuse following current in AR 10988 
and no significant current structure in AR 10989. These maps also show that 



^http: / / www.lmsal.com/solarsoft /latest_events 
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Table 1. Properties of the field measurements and NLFFF models of ARs 10987-9 in 
Figures riH13l The parameter CWsin is the current-weighted average of the sine of the 
angle between the current density and the magnetic field. 



Total magnetic fiux Ftot (Mx) 
Net magnetic fiux Fnet (Mx) 
Potential field magnetic energy Wj^ (erg) 
NLFFF magnetic energy VK[B] (erg) 
Normalized magnetic energy VK[B]/VK7r 
Relative magnetic helicity AH (G^cm^) 
Characteristic helicity AH/F^^^ 
CWsin (see caption) 



AR 10987 


AR 10988 


AR 10989 


2.3 X 10^2 


2.2 X 1022 


1.1 X 1022 


-8.0 X 10^1 


4.1 X 1021 


9.7 X 1020 


2.9 X 10^2 


3.8 X 10^2 


6.0 X 10^1 


3.1 X 10^2 


4.0 X 10^2 


6.1 X 10^1 


1.07 


1.05 


1.02 


2.9 X 10^2 


2.0 X 10^2 


2.2 X 10^1 


5.5 X 10-3 


4.1 X 10-3 


1.7 X 10-3 


0.05 


0.07 


0.1 



the vertical electric current is well balanced in each case. Therefore the electrical 
currents must flow in both directions along the field trajectories, from positive to 
negative polarity and vice versa, in approximately equal quantities. Connection 
due to flux imbalance thus corresponds to non current carrying fleld lines between 
the three active regions. 

Figures [TTIfTSl shows plots of representative fleld lines from nonlinear force- free 
models of the three regions. These plots contain the same fleld of view as the plots 
in Figures [SlfTOl The model for AR 10987 includes some twisted lines connecting 
the two sunspots while the models for AR 10988 and AR 10989 appear less 
twisted. Like AR 10987, AR 10988 has a simple bipolar appearance with long 
loops connecting fairly localized polarities. AR 10989, on the other hand, has 
more dispersed magnetic polarities and many short loops distributed over a large 
area. Table [T] shows selected physical parameter values for the magnetograms 
and models. The magnetic energy for a fleld B is deflned by the volume integral 
J (|Bp)/87r(iy. and W [B] are the energies associated with the unique poten- 
tial fleld Bt^ determined from the boundary data for the vertical fleld component 
and the NLFF fleld B, and the difference between them is the free magnetic 
energy, i.e., the maximum energy that can be released during a flare and/or 
a CME. ARs 10987 and 10988 have approximately equal unsigned magnetic 
flux, and signiflcantly more than AR 10989. On the other hand, the model for 
AR 10988 has the most magnetic energy, signiflcantly more than the model for 
AR 10987, and the model for AR 10989 has the least. We estimate from the PFSS 
model for CR 2068 that the global coronal fleld (out to 2.5 solar radii) contains 
about 1.8 X 10^^ ergs of magnetic energy, about 40% of which is accounted for 
by the NLFFF models for the active regions. AR 10989 also has signiflcantly 
less free magnetic energy than the other two regions which have approximately 
equal free magnetic energy. AR 10987 does, however, contain proportionally the 
largest free magnetic energy. AR 10987 also contains signiflcantly more relative 
magnetic helicity than AR 10988 and AR 10989 contains very little relative 
magnetic helicity. The relative magnetic helicity can easily be computed using 
XTRAPOL, which solves for the vector potential A, via the volume integral 
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Ai7 = J{A^A^)-{B-B^)dV ( [Finn and Antonsen, 1985D . The characteristic 
hehcity ( [Welsch and McTiernan, 2010 ), the relative magnetic hehcity normal- 
ized by the square of the mean magnetic flux, shows this difference between the 
regions more starkly. The field of the AR 10987 model does indeed appear to be 
more twisted than the fields of the other two models. Figure [H] shows a subset 
of field lines whose foot-points are located in a strong current concentration in 
the negative sunspot. These twisted, S-shaped lines highlight the twist in this 
active region's field structure. According to the NLFFF models the other two 
regions did not include such twisted structure. Nevertheless, according to the 
models, none of these regions contains enough free magnetic energy to power a 
major flare. 



Wheatland, Sturrock, and Roumeliotis (j2000p introduced a parameter mea- 



suring the force-freeness of a numerical model, also included in Table [H The 
quantity CWsin is the current- weighted average of the sine of the angle between 
the current density and the magnetic field. This parameter is equal to zero for 
an exact force- free magnetic field. For the three models for ARs 10987-9 plotted 
in Figures [TTlfT3l the CWsin parameter took reasonable values. 

The Hinode XRT soft X-ray image of AR 10987 on March 27 shown in Fig- 
ure [15] features an S-shaped structure commonly referred to as a sigmoid. The 
magnetic field model plots in Figure [TTI show evidence of a twisted field joining 
the two sunspots, with twist of the same handedness as the structure in the 
Hinode image. However, the plots of the modeled field also shows much structure 
that does not appear in the Hinode image. If one plots only field lines with 
foot points located within an intense concentration of vertical current, shown in 
Figure [HI an S-shaped structure stands out, bearing a striking resemblance to 
the S-shaped structure in the Hinode image. The X-ray image appears brightest 
where the twist appears greatest and where the overlying field lines meet the 
twisted structure at a large angle. 

The Hinode XRT image of AR 10988 on March 29 in Figure [15] shows plasma 
structure fanning out from the leading sunspot to the more diffuse following flux. 
In soft X-rays AR 10987 appears more compact than AR 10988 and AR 10989. 
In the plots of the modeled field for AR 10987 in Figure [121 the asymmetric 
structure of the short field lines reproduces features of the Hinode image. For 
example, the lines ending at the South edge of the leading sunspot originate from 
locations due South of the spot, corresponding to a distinct lobe at the South- 
West of the plasma structure seen in the Hinode image. Meanwhile, structures 
ending at the North, East and South-East edges of the spot curve in from the 
East in a manner common to the Hinode image and the model. 

Despite the lack of significant current structure in AR 10989, the Hinode XRT 
image of AR 10989 on March 31 shows contiguous bright structure. Its brightness 
distribution in soft X-rays is more compact than that of AR 10988 but less so 
than that of AR 10987. The set of short closed field lines connecting intense field 
concentrations in Figure [13] together resembles the bright X-ray structure. 

The XRT is sensitive to a broad range of temperatures, peaking in sensitivity 
between about 6-8 MK. For further information we also plot in Figures [T6l[T8] 
images from STEREO/SECCHI/EUVI with three narrow-band filters, 171 A 
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Figure 15. Hinode XRT images of AR 10987 on 2008 March 27 at 0617 UT (top left), 
AR 10988 on 2008 March 29 at 1759 UT (top right) and AR 10989 on 2008 March 31 at 
1802 UT (bottom). 



(corresponding to about 1 MK), 195 A (about 1.5-1.6 MK) and 284 A (about 
2.5 MK). 

Figure [m shows in the top pictures the STEREO/SECCHI/EUVI 171 A 
ahead and behind images of AR 10987 on March 27 and Figure [19] shows the 
vertical integral of the electric current squared J • J. There is no sign of twisted 
S-shaped structure in any of these plots. It appears that the sigmoidal structure 
in the Hinode XRT image is not well populated with plasma emitting at EUV 
wavelengths. The middle pictures of Figure [16] shows the STEREO behind and 
ahead images of AR 10988 on March 29. There is some resemblance between the 
modeled field plots in Figure [H and STEREO images of AR 10988. While there 
are clearly significant misalignment angles between the modeled field trajectories 
and the plasma loops there are striking qualitative similarities between them. 
This might suggest that most of the AR 10987 field is populated with X-ray 
plasma whereas AR 10988 contains plasma at a greater range of EUV and X- 
ray temperatures. Finally, the bottom pictures of Figure [TH the STEREO EUV 
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Figure 16. STEREO/SECCHI/EUVI 171 A images for AR 10987 on 2008 March 27 at 
1536 UT (top pictures), AR 10988 on 2008 March 29 at 1528 and 1526 UT (middle pictures) 
and AR 10989 on 2008 March 31 at 1741 UT (bottom pictures) from the ahead (left pictures) 
and behind (right pictures) spacecraft. 
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Figure 17. STEREO/SECCHI/EUVI 195 A images for AR 10987 on 2008 March 27 at 
1535 UT (top pictures), AR 10988 on 2008 March 29 at 1521 UT (middle pictures) and 
AR 10989 on 2008 March 31 at 1745 UT (bottom pictures) from the ahead (left pictures) and 
behind (right pictures) spacecraft. 
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Figure 18. STEREO/SECCHI/EUVI 284 A images for AR 10987 on 2008 March 27 at 
1546 UT (top pictures), AR 10988 on 2008 March 29 at 1526 UT (middle pictures) and 
AR 10989 on 2008 March 31 at 1746 UT (bottom pictures) from the ahead (left pictures) and 
behind (right pictures) spacecraft. 
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171 A images of AR 10989, do not reveal so much organized loop structure in 
AR 10989 as in AR 10987 and AR 10988. 

Figure [m shows the corresponding STEREO/SECCHI/EUVI 195 A images. 
The main difference is that ARs 10987 and 10989 appear brighter in this tem- 
perature range. More of the peripheral structure of these regions is visible at 
195 A than at 171 A , and both are brighter than AR 10988. The 195 A images 
of AR 10988 appear very similar to the 171 A images, retaining the resemblance 
to the field structure shown in Figure [T2j There is a hint of sigmoidal structure 
in the behind image of AR 10987 but it is not a clear sigmoidal signature. 

The near- simultaneous STEREO/SECCHI/EUVI 284 A images are plotted 
in Figure [181 These images differ significantly from those in Figures [16] and 
[TTl Here the structure is more diffuse and only the structure at the cores of 
the active regions are generally visible. In fact these images are more like the 
soft X-ray images from Hinode shown in Figure [TH At 284 A AR 10987 has a 
pronounced sigmoidal structure very similar to appearance in soft X-rays. The 
asymmetric loop system of AR 10988 is clearly delineated in the 284 A behind 
plot, the middle right plot of Figure [181 ^ind in the Hinode image, the top right 

plot of [m 

According to the Hinode XRT and STEREO EUVI images, then, AR 10987 
has a sigmoidal structure that is populated with hot plasma and this structure 
corresponds to a set of 'S '-shaped current-carrying field lines in the NLFFF 
model. The surrounding fields in the model shows little resemblance to the 
structure in the cooler EUVI images. AR 10988 has an asymmetric loop system 
that is visible at all temperatures and structures leaving the field of view to the 
East and West visible at the cooler temperatures. These structures correspond 
reasonably well to features in the NLFFF model. AR 10989 is much less coher- 
ently structured than the other two regions but is visible at all temperatures, 
although its peripheral structures are only visible at the cooler temperatures. Its 
basic shape resembles the set of the most intense field structures of the model 
but there is little correspondence between individual plasma loops and modeled 
field trajectories. The other two regions had significantly more intense fields than 
this one. 

Figure [19] shows the vertical integral of the modeled coronal current density 
squared J • J for ARs 10987 and 10988 both on March 27. The vertically inte- 
grated current density squared (or ohmic heating for a uniform resistor) does not 
forward-model actual soft X-ray or EUV emission but it is often used as a proxy 
for a qualitative comparison. Comparing the AR 10987 plots to Figures [TTl E] 
and [15] the integrated current structure shows no sign of the twisted, sigmoidal 
structure of the field trajectories and XRT brightness pattern. Indeed, none of 
the models for AR 10987 show sigmoidal structure in their integrated vertical 
current profiles. The twisted current-carrying fields seen in Figure [TH perhaps 
extend to high into the atmosphere and are perhaps too weak to compete with 
the lower- lying strong field that dominates the integrated current image. The plot 
for AR 10988, however, does resemble the structure of the field in Figures [TTl 
and the EUV emission patterns in the middle pictures of Figures [16] and [17] On 
the other days during its disk passage the current structure of AR 10988 looks 
similar but some of the details are missing. 
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Figure 19. The vertically integrated modeled coronal current densities for ARs 10987 (top) 
and 10988 (bottom) on March 27. 
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Figure 20. STEREO/SECCHI/EUVI 171 A image of the three active regions 10987-9 on 
2008 March 29 at 0409 UT from the behind spacecraft. 

The STEREO/EUVI 171 A image in Figure [20] shows evidence that the three 
active regions are magneticahy connected. In particular, AR 10989, which is 
approximately flux-balanced, has structures extending outwards both to the East 
and to the West. 

As the top panel of Figure [211 shows, the PFSS model for the global coronal 
field of rotation 2068 also indicates that all three active regions are magneti- 
cally connected. In the model, magnetic fiux directed towards the West side of 
AR 10989 originates from the East side of AR 10988 and, the fiux of AR 10989 
being nearly balanced, positive fiux from the East side of AR 10989 connects 
to quiet-Sun fields to the East. The connectivities of AR 10987 and AR 10988 
are more unusual. Because these regions have significant positive and negative 
fiux imbalances, respectively, there is less connection between their neighboring 
positive and negative polarities than between their more separated negative and 
positive fields. Flux connecting these fields must loop over both active regions. 
The second plot shows the same model but also includes coronal holes and 
streamer-belt field lines. Comparing the two plots, the tall, looping field lines 
connecting AR 10987 and AR 10988 are contained within the streamer belt while 
a significant quantity of AR 10987's fiux, directed towards the Western part of 
the region, is open field. All but this Western part of AR 10987 is contained 
within the streamer belt. Almost all of AR 10988 is also within the streamer 
belt but the part outside, the Eastern edge, is mostly connected to AR 10989 
to the East. AR 10989 is unique in lying entirely outside albeit close to the 
streamer belt. An active region, consisting of closed loop structures, located 
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Figure 21. Two pictures of the PFSS model for CR 2068. The two pictures show the synoptic 
magnetogram in greyscale. The top picture shows field lines with foot-point fields of strength 
15 G or higher in yellow. The bottom picture shows field lines with foot-point fields of strength 
45 G or higher in yellow, streamer-belt field lines in blue and positive and negative coronal 
holes in red and green dots. 



mostly outside the streamer belt among ambient open field may be expected to 
have more fields in contact with topologically different domains and therefore 
be more flare-productive than a region embedded within closed fields under 
the streamer belt. This may help to explain why, although this region has the 
least free magnetic energy and helicity of the three regions, it is also the most 
flare-productive ( |Webb et al.^ 20101 ). 

Figure [22] shows the net, positive and negative magnetic fluxes of the three 
active regions over five daily observations. Almost all measurements show a sig- 
nificant positive flux imbalance in AR 10987, a significant negative flux imbalance 
in AR 10988 and approximately balanced flux in AR 10989. The measurements 
for AR 10987 seem to be approximately steady in time whereas AR 10988 and 
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Figure 22. Temporal evolution over five-day intervals of the magnetic (left plots) and electric 
current (right plots) fluxes for AR 10987 (top plots) 10988 (middle plots) and 10989 (bottom 
plots). Red/blue/black symbols represent positive/negative/net fluxes. 
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Figure 23. Temporal evolution over five-day intervals of the modeled magnetic energies (left 
plots) and relative magnetic helicities (right plots) for AR 10987 (top plots) 10988 (middle 
plots) and 10989 (bottom plots). Crosses represent NLFFF models and circles represent the 
corresponding potential-field (current-free) models. The potential-field models have zero rela- 
tive magnetic helicity. Red/blue symbols represent models integrated from positive/ negative 
vertical magnetic fields. 
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AR 10989 show significant variation. While AR 10988's negative flux declines 
steadily in the plots, its positive flux has a large peak on March 30th, resulting in 
a particularly large flux imbalance on that day. Inspection of the images reveals 
anomalously strong plage field measurements on that day that seem unlikely to 
be accurate. The positive and negative fluxes of AR 10989 increase in size with 
distance from disk-center. This apparent dependence on position on the solar 
disk is most likely an artifact due to noise. The further the region is displaced 
from disk-center, the more the signal for Bz derives from transverse as opposed 
to longitudinal measurements. Transverse field measurements are intrinsically 
noisier than longitudinal field measurements. It seems that the measurements 
for AR 10989 are more affected by this noise than are the measurements for the 
other two regions because AR 10989's fields are relatively weak. 

Also plotted are the net, positive and negative electric current fluxes of 
the regions. The current fluxes are very well balanced in these measurements. 
AR 10987 has a gradual overall increase in current. Moreover, the current is 
well balanced at each sunspot, suggesting that local twisting motions are adding 
twist to the field. AR 10988's graphs are more ambiguous. The electric current 
flux of AR 10989 shows the same dependence on position on the solar disk as 
its magnetic flux. 

Figure [23] shows plots of the free magnetic energies and helicities of the models 
of the three regions based on the same five-day daily observations. The plots 
include models based on both positive- and negative-field boundary conditions. 
As we mentioned in Section [2] the reconstruction method can use boundary 
conditions for a either where the vertical field is positive or where it is negative. 
While innovative methods have been developed to take into account both po- 
larities by [Wheatland and Regnier ()2009p and Amari and Aly[ (j2010p , our active 



regions are nearly potential so that the positive- and negative-polarity solutions 
tend to match quite well. When they do not match, there is generally a simple 
explanation. Hence for simplicity we retain the simpler approach of integrating 
from each polarity separately. In Figure [23] the left plots of magnetic energy 
also include potential-field models for comparison. It is clear that none of these 
models is far from a potential state. AR 10988 has the most photospheric flux of 
the three regions and also has the most magnetic energy. All three regions show 
decreasing trends in magnetic energy, with outlying data points for AR 10988 
on March 30 and for AR 10989 on March 31 interrupting these trends. All 
three regions appear to have decreasing free magnetic energy, perhaps due to 
cumulative effects of turbulent diffusion. It is notable that the model parameters 
for AR 10989 do not show the same dependence on position on the solar disk 
as the magnetic and electric current fluxes do for that region. AR 10987 has 
proportionally significantly more free magnetic energy than the other regions, 
reflected in its twisted, sigmoidal structure. However, the free energy contained 
in even this region is much less than would be needed to power a major flare. 

The graphs of magnetic helicity are more complicated because this quantity 
is very sensitive to problems in the data and model. In extreme cases the model 
can even give the wrong sign for the helicity. For example, the positive-polarity 
model for AR 10988 on March 27 had a problem because the large positive 
flux balance of the region caused a significant amount of current-carrying flux 
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to leave the numerical domain, instead of closing within the region. Because 
magnetic helicity depends on the twist and interlinking of fields, such a prob- 
lem can significantly affect the relative helicity of the modeled field. Similar 
problems occur with the negative-polarity model for AR 10987 on March 25th 
and 28th and for the positive-polarity model for AR 10988 on March 29th. In 
such cases it is clear that the anomalous models should be disregarded. Judging 
from those models that are free of such problems, whose helicities are plotted 
in the right half of Figure [23l ARs 10987 and 10988 most likely have positive 
magnetic helicity while AR 10989 does not have significant helicity. These results 
are consistent with the statistical result ([Pevtsov, Canfield, and Metcalf, 1995| 



Pevtsov, Canfield, and Latushko, 2001D that active regions in the northern/southern| 



hemisphere tend to have negative/positive helicity. 

IWelsch and McTiernan (2010] calculated from time series of MDI longitudinal 
magnetograms that AR 10988 contains the most magnetic fiux and AR 10989 
the least. They find that, while the three regions are fed similar quantities of 
relative helicity fiux, the characteristic magnetic helicity fiux is greatest for 
AR 10989 because of its relatively small magnetic fiux. This is the only evidence 
they derived from the MDI magnetograms that AR 10989 should have been the 
most active of the regions. As Welsch and McTiernanl(|2Q10j) have noted, we find 
little evidence of significant relative magnetic helicity in the NLFFF models for 
AR 10989. We already know that AR 10989 was not a very structured region and 
that it lay mostly in a region of open field outside the streamer belt. Perhaps 
the relative magnetic helicity transported into the corona was not effectively 
trapped in the coronal field of AR 10989. 



5. Conclusion 

We have calculated PFSS models for Carrington rotations 2067-9 and NLFFF 
models for the three active regions of the WHI, NOAA ARs 10987-9. Together 
the models reveal the following properties of coronal magnetic field structure 
around the time of the WHI. 

The global coronal field structure was very asymmetric around the time of 
the WHI with all significant activity crowded between about 180° and 270° of 
Carrington longitude. This lopsided activity distribution, combined with rela- 
tively weak polar fields, resulted in a global coronal field structure that was 
far from the nearly axisymmetric, dipolar configuration expected during solar 
minimum. There was a large warp in the streamer belt centered at around 240° 
of Carrington longitude and a displacement of the streamer belt from the equator 
of 30 to 40° during all three rotations CR 2067-9 in both the modeled fields and 
the STEREO/SECCHI/CORl streamer observations. Furthermore, much of the 
fiux connecting to the ecliptic plane originated from low-latitude coronal holes, 
an unusual state of affairs for a near-minimum corona. 

During the WHI itself, CR 2068, the magnetic activity appeared to be par- 
ticularly well organized in three active regions about 30° apart in longitude and 
approximately equal latitude of about —10°. All three were bipolar and had 
leading negative polarity, consistent with Hale's law, although their tilt angles 
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appeared to be different. There were further significant differences between the 
ARs. The photospheric field of AR 10987 was dominated by two sunspots of 
approximately equal strength. AR 10988 had a strong leading sunspot but widely 
distributed, less intense following flux while AR 10989 had no sunspot and was 
significantly weaker and more fragmented than the other two. ARs 10987 and 
10988 had comparable magnetic fluxes and AR 10989 had significantly less. 
According to NLFFF models AR 10988 had the largest magnetic energy of 
the three regions although AR 10987 had approximately the same amount of 
photospheric magnetic flux and had proportionally more free magnetic energy. 
AR 10987 also had the most relative magnetic helicity of the three regions 
and AR 10989 the least. The models for ARs 10987 and 10988 indicate posi- 
tive relative magnetic helicity, consistent with a well-known statistical tendency 
for southern- hemisphere active regions, while the relative magnetic helicity of 
AR 10989 was not significant and its sign was not determined by the models. 
None of the regions had free magnetic energy or helicity normally associated 
with a flaring active region. However, AR 10987 had a coherent S-shaped or 
sigmoidal structure connecting its two sunspots. This was embedded in a bipolar 
loop system in the NLFFF model and had foot points in an intense current 
concentration in the positive-polarity sunspot. This structure was also evident 
in the high-temperature STEREO/SECCHI/EUVI and Hinode/XRT images. 
ARs 10988 and 10989 had no corresponding structure in the models or the 
observations. The PFSS model for CR 2068 showed that the three regions were 
highly interconnected, and that ARs 10987 and 10988 were located mostly within 
the helmet streamer belt whereas most of AR 10989 lay outside the streamer 
belt among open fields. 

According to Webb et aU ([2QTq1) AR 10989 was the most CME- and flare- 
productive of the three regions. While IWelsch and McTiernanl(|2010p argue that 
the flare activity may have been more balanced among the three regions, sev- 
eral predictors of flare activity that they calculated based on the photospheric 
magnetic field indicated that AR 10988 should have been the most active re- 
gion of the three and AR 10989 the least. Of the magnetic parameters they 
studied, the characteristic helicity, the cumulative relative magnetic helicity flux 
normalized by the square of the mean magnetic flux, was the only one for which 
AR 10989 had a higher value than the other two regions. Our NLFFF models 
tell a similar story. Of the three regions, AR 10989 had the lowest free magnetic 
energy and relative magnetic helicity. In the present paper we have computed 
the two solutions each one corresponding to injection of electric current from one 
polarity. It would be interesting to extend this study using the new Optimization 
Grad Rubin Method (OGRM) described in Amari and Aly (2010) for the NLFF 
calculation to find a unique solution obtained by taking into account electric 
current from both polarities. 

IWebb et aL\ (f2QTQ|) and IWelsch and McTiernanl (|2QTq1) discuss possible rea- 
sons why AR 10989 was relatively active. According to helioseismological calcu- 
lations, subsurface vortical flows were occurring beneath all three active regions 
but AR 10989 was distinguished in having strong neighboring vortical flows of op- 



posite sign, a characteristic of flare-productive active regions (Webb et a/., 2010). 
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Moreover the field of AR 10989 may have evolved between the time of the mag- 
netic observations near disk-center and the time of the CME activity observed 
near AR 10989 at the limbs (Welsch and McTiernan, 2010D. Our PFSS model in- 



dicates that AR 10989 was situated among open fields while our NLFFF models 
show how relatively weak and unstructured the field of AR 10989 was compared 
to the other two regions. One tentative explanation is that the equilibrium of 
AR 10989's field was relatively unstable. Perhaps the subsurface vortical motions 
and photospheric helicity flux easily destabilized the coronal field, whereupon 
magnetic energy and helicity readily escaped from the vicinity of AR 10989 via 
the ambient open fields. What is more certain is that we have much to learn 
about the interacting processes of solar activity. 
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